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Abstract 

The incidence of Dengue epidemiologic disease has grown in recent decades. In this paper 
an application of optimal control in Dengue epidemics is presented. The mathematical model 
includes the dynamic of Dengue mosquito, the affected persons, the people's motivation to 
combat the mosquito and the inherent social cost of the disease, such as cost with ill indi- 
viduals, educations and sanitary campaigns. The dynamic model presents a set of nonlinear 
ordinary differential equations. The problem was discretized through Euler and Runge Kutta 
schemes, and solved using nonlinear optimization packages. The computational results as well 
as the main conclusions are shown. 
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1 Introduction 

Dengue is a mosquito mostly found in tropical and sub-tropical climates worldwide, mostly in 
urban and semi-urban areas. It can provoque a mosquito-borne infection, that causes a severe 
flu-like illness, and sometimes a potentially lethal complication called dengue haemorrhagic fever, 
and about 40% of the world's population are now at risk. 

The aim of this paper is to present an attempt to apply quantitative methods in the opti- 
mization of investments in the control of Epidemiologic diseases, in order to obtain a maximum 
of benefits from a fixed amount of financial resources. This model includes the dynamic of the 
growing of the mosquito, but also the efforts of the public management to motivate the population 
to break the reproduction cycle of the mosquitoes by avoiding the accumulation of still water in 
open-air recipients and spraying potential zones of reproduction. 

The paper is organized as follows. Next section presents the dynamic model for dengue epi- 
demics, where the variables, parameters and the control system are defined. Then, the numerical 
implementation and the strategies used to solve the problem are reported. Finally, the numerical 
results are presented and some conclusions are taken. 



*This is a preprint of a paper whose final and definite form has been published in American Institute of Physics 
Conf. Proc. 1168 (2009), no. 1, 1385-1388. DDI : 10. 1063/1 .3241345 
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2 Dynamic model 



The model described in this paper is based on the model proposed in [3]. 
The notation used in the mathematical model is as follows. 

State Variables: 

xi {t) density of mosquitoes; 

X2{t) density of mosquitoes carrying the virus; 

X3{t) number of persons with the disease; 

X4{t) level of popular motivation to combat mosquitoes (goodwill). 

Control Variables: 

ui{t) investments in insecticides; 

U2{t) investments in educational campaigns. 

Parameters: 

an average reproduction rate of mosquitoes; 
aM mortality rate of mosquitoes; 

/? probability of contact between non-carrier mosquitoes and affected persons; 
77 rate of treatment of affected persons; 

H amplitude of seasonal oscillation in the reproduction rate of mosquitoes; 
p probability of persons becoming infected; 

9 fear factor, reflecting the increase in the population's willingness to take actions 
to combat the mosquitoes as a consequence of the high prevalence of the disease 
in the specific social environment; 

T forgetting rate for goodwill of the target population; 

(f phase angle to adjust the peak season for mosquitoes; 

w angular frequency of the mosquitoes's proliferation cycle, corresponding to 52 
weeks period; 

P population in the risk area (usually normalized to yield P = 1); 
7£) the instantaneous costs due to the existence of affected persons; 
7i? the costs of each operation of spraying insecticides; 
the cost associated to the instructive campaigns. 

The Dengue epidemic can be modeled by the following nonlinear time- varying state equations. 
Equation ([!]) represents the variation of the density of mosquitoes per unit time to the natural 
cycle of reproduction and mortality (a/j and q;m), due to seasonal effects fisin(ujt + ip) and to 
human interference —Xi{t) and ui{t): 

(ix\ 

— = [an (1 - nsin{ujt + ip)) ~ um - X4.{t)] xi{t) - ui{t). (1) 

Equation ^ expresses the variation of the density X2 of mosquitoes carrying the virus. The 
term [an (1 — iJ,sin{ijjt + tp)) — aM — 2^4 (i)] X2(t) represents the rate of the infected mosquitoes and 
f3[xi{t) —X2{t)]x3{t) represents the increase rate of the infected mosquitoes due to the possible 
contact between the non infected mosquitoes xi(t) — X2(t) and the number persons with disease 
denoted by X3: 



dx2 

—— = [ofK (1 - ^isin{ujt + ip)) - aM ~ x^it)] X2{t) + I3[xi{t) - X2(t)] X3{t) - ui(t). (2) 
at 

The dynamics of the infectious transmission is presented in equation fS]). The term —riX3{t) 
represents the rate of cure and px2{t) [P — X3{t)\ represents the rate at which new cases spring up. 
The factor [P — X'i{t)\ is the number of persons in the area, that are not infected: 

^ = -77X3 {t) + pX2 (t) [P - X3 (t)] . (3) 
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Equation (U) is a model for the level of popular motivation (or goodwill) to combat the re- 
productive cycle of mosquitoes. Along the time, the level of people's motivation changes and, as 
consequence, it is necessary to invest in educational campaigns designed to increase consciousness 
of the population under risk by a proper understanding of the determinants involved with specific 
disease. The expression —TX4{t) represents the decay of the people's motivation with time, due to 
forgetting. The expression 6x3 (t) represents the natural sensibilities of the public due to increase 
in the prevalence of the disease. 

^ = -TX^it) + eX3{t) + U2{t). (4) 

The goal of the problem is to minimize the cost functional 

J[ui{.),U2{.)]^ {jDxlit) + jFulit) + jEulit)}dt. (5) 

Jo 

This functional includes the social costs related to the existence of ill persons, jax^it), the 
recourses needed for spraying of insecticides operations, jpuKt), and for educational campaigns, 
7£;M|(t). The model for the social cost is based in the concept of goodwill explored by Nerlove 
and Arrow '5'. 

Due to computational issues, the optimal control problem (IT])-(IS|) that is in the Lagrange form, 
was converted into an equivalent Mayer form. Hence, using a standard procedure to rewrite the 
cost functional [6 , the state vector was augmented by an extra component X5, 

^ = Jd4 it) + 1FU\ (t) + ^EUI it) (6) 

leading to the equivalent terminal cost problem of minimizing 
with given tf, subject to the control system and (jH]). 



3 Numerical implementation 

The simulations were carried out using the following normalized numerical values: a^. — 0.20, 
an = 0.18, /3 = 0.3, 77 = 0.15, fi = 0.1, p = 0.1, 9 = 0.05, r = 0.1, = 0, w = 27r/52, P = 1.0, 

= 1.0, -fF = 0.4, 7B = 0.8, xiiO) = 1.0, X2i0) = 0.12, xs^O) = 0.004, and 2:4(0) = 0.05. These 
values are available in the paper |3^. The final time used was i/ = 52 weeks. 

To solve this problem it was necessary to discretize the problem. Two methods were selected: 
a first order, the Euler's scheme, and a Runge Kutta's sheme of second order [1]. In both cases, 
it is assumed that the time t = nh moves ahead in uniform steps of length h. If a differential 
dx 

equation is written like — = f(t, x), it is possible to make a convenient approximation of this. In 
dt 

the Euler's scheme the update is given by 
while in the Runge Kutta's method is 

Xn+l ~ + I [f{tn,Xn) + f Xn+l)] ■ 

This approximation Xn+i of x{t) at the point has an error depending on and h^, 

for the Euler and Runge Kutta methods, respectively. This discretization process transforms 
the dengue epidemics problem into a standard nonlinear optimization problem (NLP), with an 
objective function and a set of nonlinear constraints. This NLP problem was codified, for both 
discretization schemes, in the AMPL modelling language [3]. 
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Euler's method Rungc Kutta's method 



Knitro 


h 


# var. 


# const. 


# iter. 


time (sec.) 




h 


# var. 


# const. 


# iter. 


time (sec.) 




0.5 
0.25 
0.125 


727 
1455 
2911 


519 
1039 
2079 


113 
68 
85 


2.090 
2.210 
7.240 


0.5 
0.25 
0.125 


728 
1456 
2912 


520 
1040 
2080 


64 
82 
70 


1.980 
5.550 
9.740 


Snopt 


h 


# var. 


# const. 


# iter. 


time (sec.) 




h 


# var. 


# const. 


# iter. 


time (sec.) 




0.5 
0.25 
0.125 


727 
1455 
2911 


519 
1039 
2079 


175 
253 
252 


4.07 
19.2 
105.4 


0.5 
0.25 
0.125 


728 
1456 
2912 


520 
1040 
2080 


223 
219 
420 


10.52 
39.7 
406.67 



Table 1: Numerical results 



Two nonlinear solvers with distinct features were selected to solve the NLP problem: the Kni- 
tro and the Snopt. The first one [2] is a software package for solving large scale mathematical 
optimization problems based mainly on the Interior Point (IP) method. Snopt [5] uses the SQP 
(Sequential Quadratic Programming) philosophy, with an augmented Lagrangian approach com- 
bining a trust region approach adapted to handle the bound constraints. The NEOS Server [7] 
platform was used as interface with both solvers. 

4 Computational results 

Table [1] reports the results for both solvers, for each discretization method using three different 
discretization steps {h = 0.5, 0.25, 0.125), rising twelve numerical experiences. The columns # var. 
and # const, mean de number of variables and constraints, respectively. The next columns refer to 
the performance measures - number of iterations and total CPU time in seconds (time for solving 
the problem, for evaluate the objective and the constraints functions and for input /output). The 
computational experiences were made in the NEOS server platform - in this way the selected 
machine to run the program remain unknown as well as its technical specifications. 

The optimal value reached was ~ 3E — 03 for all tests. Comparing the general behaviour of the 
solvers one can conclude that the IP based method (Knitro) presents much better performance than 
the SQP method (Snopt) in terms of the measures used. Regarding the Knitro results, one realize 
that the Euler's discretization scheme has better times for h = 0.25 and h = 0.125 and similar time 
for h = 0.5, when compared to Runge-Kutta's method. Another obvious finding, for both solvers, 
is that the CPU time increases as far as the problem dimension increases (number of variables 
and constraints). With respect to the number of iterations, Snopt presents more iterations as 
the problem dimension increases. However this conclusion cannot be taken for Knitro - in fact, 
doesn't exist a relation between the problem dimension and the number of iterations. The best 
version tested was Knitro using Runge-Kutta with h = 0.5 (best CPU time and fewer iterations), 
and the second one was Knitro with Euler's method using h = 0.25. An important evidence of this 
numerical experience is that it is not worth the reduction of the discretization step size because 
no significative advantages are obtained. 

5 Conclusions 

We solved successfully an optimal control problem by direct methods using nonlinear optimiza- 
tion software based on IP and SQP approaches. The effort of the implementation of higher order 
discretization methods brings no advantages. The reduction of the discretization step and conse- 
quently the increase of the number of variables and constraints doesn't improve the performance 
with respect to the CPU time and to the number of iterations. We can point out the robustness 
of both solvers in spite of the dimension problem increase. 
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